Direct measurements of the dynamical correlation length indicate its divergence at an 

athermal glass transition 
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The super-cooled iV3 model exhibits an increasingly slow dynamics as density approaches the 
model's random closest packing density. Here, we present a direct measurement of the dynamical 
correlation function Gt{r,t), showing the emergence of a growing length scale £4 across which the 
dynamics is correlated. The correlation length measured, up to 120 lattice sites, power-law diverges 
as the density approaches p t , the density at which the fluid phase of the model is predicted to 
terminate. It is shown that the four-point susceptibility, often used as an agent to estimate £4, does 
not depend simply on the latter. Rather, it depends strongly on the short-range behavior of Gi{r, t). 
Consequently, \i peaks before £4 reaches its maximal value. The two quantities should therefore be 
studied independently. 



Dynamics of super-cooled liquid glass-formers is char- 
acterized by non-exponential temporal relaxation (see re- 
cent review [l|). Two natural, but fundamentally dif- 
ferent, explanations can be put forward. Relaxation 
might be locally exponential, but the typical relaxation 
timescale varies spatially. Global response functions be- 
come non-exponential upon spatial averaging, due to 
the spatial distribution of relaxation times. Alterna- 
tively, relaxation might be complicated and inherently 
non-exponential, even locally. Experimental and theo- 
retical works 0] suggest that while the latter mechanism 
is likely to contribute, relaxation (at times of order of the 
global relaxation time or shorter) is indeed spatially het- 
erogeneous. At an ideal glass transition, the relaxation 
time diverges and the dynamics remains heterogeneous 
for arbitrary long times. The study of dynamical het- 
erogeneities has attracted much interest recently 
as part of the attempt to decipher the mystery of glassy 
dynamics. 



A physical characterization of dynamic heterogeneity 
entails the determination of the typical lifetime of the 
heterogeneities, as well as their typical lengthscale. A 
clear demonstration of the heterogeneous dynamics is ob- 
tained by tracking the mobile particles at a given time. 
One observes clusters of highly mobile particles as well 
as clusters of particles barely moving at all [IM2- The 
typical size of these clusters is quantified by the spatial 
decay of the mobility correlation function. One there- 
fore looks for the correlation between the displacements 
over a time interval t of particles at mutual distance r. 
This mobility-mobility correlation function was first in- 
troduced in [13], as a tool to discover cooperative re- 
gions in numerical simulations of glass-forming liquids. 
However, the first attempts to analyze these correlation 
function [I3} were very limited, and thus showed no in- 
creasing correlation length. It turns out that the essential 
ingredient to find the dynamical correlations is not really 
the mobility, but rather the fact that we are calculating a 
four-point correlation function, in contrast with standard 



two-point functions. The four point density correlation 



G±(r,t) = (p(0,Q)p(0,t)p(r,0)p(r,t)) 

-(p(0,0)p(0,t))(p(r,0)p(r,t)), 



(1) 



was first suggested in [14j in order to look for a growing 
correlation length. Yet, direct measurements of G4 turn 
out to be technically demanding. Instead, it is common 
to study its spatial integral, which is nothing but the 
fluctuations of the two-point correlation function C (t) = 
(p(0, t)p(r, t)) — (p(0, t))(p(r, t)), or the four-susceptibility 



X4 (t) = / G 4 (r,t)dr = N[(C(t) 2 ) - (C(t)) 2 



(2) 



The four-susceptibility was indeed shown to exhibit an 
appreciable increase with the waiting time (lij ]. which 
was interpreted as a sign for an increase of the dynami- 
cal correlation length characterizing the decay of G4 (r, t) . 
The four-point susceptibility was later used extensively to 
measure dynamical heterogeneities approaching the glass 
transition and to imply the existence of growing correla- 
tion length. The time-dependence of Xi(t) was suggested 
to indicate the mechanism underlying the transition Q. 

Various theories lead to different relations between the 
correlation length £4, characterizing the spatial decay of 
G*4(r, t), and \4 UM> Thus, direct measurements 

of dynamical correlations can play a vital role in dis- 
tinguishing the various theoretical approaches. In fact, 
it was pointed out in 0, EH that X4 and £4 might de- 
pend differently on time. However the calculation of dy- 
namic correlations is hard especially in numerical simu- 
lations, because very large systems are needed in order 
to determine £4 unambiguously. Previous measurement 
of dynamic correlations were performed indirectly by ei- 
ther scaling the binder cumulant at T4 (time at which 
(t) peaks) [9J or fitting the low q behavior of Gn(q,t) 
. These measurements have shown moderate corre- 
lation lengths, up to 10 inter-particle distances. Clearly, 
these do not describe very well the behavior at a regime 
in which correlation lengths (are expected to) diverge. 
Here we report measurements of dynamical correlation 
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FIG. 1. (Color Online) Random packing of the N3 model as 
produced by the cooling protocol, density set to 0.17. The 
few particles that are movable at the end of the cooling phase 
are marked by a gray spot at their center. 



based on the evaluation of the four-point density correla- 
tion Gi{r,t) throughout the relaxation process. Having 
measured G^ir, t) directly, the correlation lengths is then 
extracted by fitting an exponentially decaying function to 
the long range behavior of Gi{r,t). 

The 7V3 model is a simple 2D model on a square lat- 
tice. Particles are cross-shaped pentamers, interacting 
only through hard-core exclusion which blocks up to the 
3 rd nearest neighbor (see Fig. [IJ. The model is known 
to undergo a first order solidification transition [19142 ll . 
where density jumps from pf ~ 0.161 to p s ~ 0.191 [21[ 
(the closest packing density is 0.2). The behavior of this 
model in the supercooled fluid regime was analyzed using 
the R matrix method 
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|23j based on the Mayer clus- 
ter coefficients of N3. The analysis predicted a critical 
termination of the supercooled fluid phase at non trivial 
density p t = 0.1717, close to the random closest pack- 
ing density of this model [24[ . Extensive MC simulations 
revealed that the dynamic relaxation of the supercooled 
fluid becomes increasingly slow as the density increases, 
diverging exactly at the density predicted by R matrix 
analysis [25j. It is therefore plausible to relate the dy- 
namical arrest in the N3 model to the predicted termi- 
nation of the thermodynamic equation of state predicted 
by the R matrix. The iV3 system is then a simple and 
convenient model-system for studies of slow dynamics in 
quenched deeply supercooled fluids. 

Here we employ the advantage of the the simple two 
dimensional N3 system and study in detail the dynamic 
correlation length near the predicted termination of the 
fluid supercooled phase. Simulations of the N3 model 
allowed for an easy and direct measurements of G4 (r, t) 
throughout the whole relaxation process for a wide range 
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FIG. 2. (Color Online) Measurements of Gi(r,t) during dy- 
namical relaxation of the N3 model, at density p — 0.1713. (a) 
t=5000 (b) t= 10000, a clear exponentially regime emerges, 
(c) t=20000, the correlation length further increases, (d) 
t=30000, further increase of the correlation length. Note that 
the amplitude of the correlation function G$(r,t) is smaller 
in this case, leading to an overall smaller integral X4 as com- 
pared to the time presented in (c) . The exponential fit in red 
may be used to extrapolate Gi{r,i) to larger distances (see 
text). 
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FIG. 3. (Color Online) Time dependence of the dynamical 
correlation length for several values of the density. The cor- 
relation length increases with time, and then saturates for as 
long as we are able to measure it. The correlation length as 
measured at T4 power-law diverges with critical exponent of 
1.0 (inset). 
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of system sizes. We performed Monte Carlo simulations 
of the model following the protocol presented in [25j |: 
simulation starts with an infinitely fast cooling, where 
particles are added whenever possible and diffuse other- 
wise, this process is stopped when the desired density is 
reached. The system is then left to relax diffusively at 
a fixed density, and measurements are performed during 
this relaxation process. Figure [2] presents the behavior 
exhibited at the different time regimes. We have studied 
lattices of sizes up to 1200 x 1200 sites in the current 
study (up to 2000 x 2000 sites in the past). Simulations 
of such large systems are crucial, as significant finite-size 
effects are shown to persist to large lattice sizes (25|. In 
concordance, we measure correlation lengths as large as 
120 lattice-sites (50 inter-particle distances). The corre- 
lation length is estimated by fitting an exponential form 
to the measured G^r.t). The only freedom in the fit is 
the spatial range upon which to fit the exponential form: 
in all of the £4 measurements below we used the range 
30 < r < 90, even though for cases with large £4 the fit 
captures the behavior far beyond this chosen range, as 
seen in Figure [2]) . 

The picture that emerges from these simulations is 
the following: at short times Gi(r,t) decays at short 
range, with no exponential tail observed. However, on 
intermediate times one clearly sees an exponential de- 
cay regime following the short-range decrease. At later 
times, (later then T4), one again sees only the short range 
behavior. Figure [3] presents the time dependence of the 
dynamic correlation length £4, for various densities of 
the N3 model. We find that the correlation length ^(t) 
grows with time, approaching a plateau. It is not possible 
to measure £4 at very long times, where G^ir, t) is very 
small beyond the short-range regime, and it is not at all 
clear whether there is an exponential regime in these long 
times. However, for as long as we are able to measure £4 
it does not seem to increase or decrease. Furthermore, 
the time dependence of xa (t) is determined by the short 
range correlation (which dictates the amplitude of the 
long range behavior) as well as the correlation length it- 
self. Due to the dependence of both £4 and the amplitude 
of Gi(r,t) on the time t, the time T4 at which Xi peaks 
differs from the time at which £4 reaches its maximal 
value. These findings highlights the non-trivial relation 
between Xi and £4, and the need to study both quantities 
separately. In addition, our measurements show that the 
dynamically correlated regions have compact (and not 
fractal) form in all times measured. 

Looking at the microscopic structure of the relaxation 
dynamics in the N3 model, one observes a clear hetero- 
geneous picture similar to that seen, for example, in the 
kinetically constrained triangular lattice gas models (1)- 
TLG and (2)-TLG 2fj. The dynamics can be described 
by a growth of mobile regions in the system, as seen in 
Figure HI until all of the particles moved when C{t) = 0. 
In course of time, the mobile regions, or clusters, grow, 



but there is no creation of new clusters. The measured 
correlation length fits this observed dynamics. The in- 
crease in the correlation length corresponds to the growth 
of clusters of mobile particles, while in later times further 
cluster growth is blocked by the neighboring clusters, re- 
sulting in a saturation of correlation length. 

Measurements of G4 have revealed an additional phe- 
nomenon. As seen in Figure [2j G^r, t) does not decay to 
at very long ranges, of order of the system's size. Close 
examination of this unexpected feature shows that the 
(unphysical) infinite-range correlation observed depends 
on system size, decreasing like 1/L 2 . It thus follows that 
its contribution to Xii t ne spatial integral over G4, per- 
sists for very large systems. It turns out that the source 
of these spurious infinite-range correlations is the initial 
configuration (that is, correlations introduced during the 
cooling protocol). To show this, we looked at the correla- 
tion function of the density p m (r, t = 0) of particles which 
are free to move at t = (end of cooling phase and start 
of the measurement time) . The global density of movable 
particles at time t = is just m = {p m {r,t — 0)), and 
the generalized susceptibility for this density is given by 

Xra = V[(m 2 ) - (TO) 2 ] 

. This quantity is closely related to X4 & t short times, 
since at these short times particles that indeed moved 
are, by and large, those that were free to move at t = 0. 
Looking at the correlation function of the movable par- 
ticles density G m (r) — (m(O)m(r)) — (m) 2 one observes 
the same infinite-range contributions. It thus follows that 
infinite-range correlations in the density of the movable 
particles emerge during the cooling protocol we used. 
Similar infinite-ranged correlations have been observed 
for other out-of-equilibrium systems [27| ■ This global ef- 
fect is then amplified in the course of dynamics, where 
movable particles grow into clusters, and affects strongly 
dynamical measurements of \4- The effect is washed out 
only for relatively long times, often of order T4 and more. 
Obviously, these infinite-range contributions to Xi have 
nothing to do with the dynamics in general, or with the 
dynamical correlation length in particular. 

Similar effects are expected in other studies, as long as 
the system is not equilibrated before the measurement 
starts. Such equilibration is always desired, but usually 
not possible when studying systems that are prone to 
crystallization (such as hard-spheres, or the N3 model). 
Thus, direct measurements of Xi m such cases are prone 
to be affected by the effect described above. However, if 
one measures Gi(r, t), a correction scheme for the initial 
condition is possible. The value of Xi> often calculated 
using N ■ Var(C(t)), does contain global contributions 
which are artifacts of the cooling protocol. Based on the 
observed dynamics, we claim that only the growth of the 
mobile regions is a local dynamical property that should 
be taken into account, and the infinite-range contribu- 
tion (present even at t = and therefore not connected 
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FIG. 4. (Color Online) Displaying relaxation dynamics of N3 
model, for p = 0.1713. Dots represent mobile particles, the 
growth of mobile regions as well as no creation of new mobile 
regions with time is observed. An area of 200x200 was chosen 
out of simulated system of 1000x1000 sites 



proached X4,max ~ {pt ~ p)~ 2 ' - The time-dependence of 
Xi{t) underscores its dependence on the short range be- 
havior of Gi(r,t) rather than the long exponential tail. 
The short range correlation determines the coefficient 
multiplying the exponential decay, and has stronger time 
dependence then and so determines Xi{t) time de- 
pendence. 

In summary, we present here a direct measurement of 
the dynamical correlation function G^{r^t), showing the 
emergence of a growing length scale £4 across which the 
dynamics is correlated. One observes a significant in- 
crease in the correlation length, up to 120 lattice sites. 
As density approaches pt, the density at which the fluid 
phase of the model is predicted to terminate, the correla- 
tion length £4 power-law diverges with (p t —p). In concor- 
dance, peak values of X4,{t) diverge, suggesting that the 
slow-down of relaxation processes near the termination 
of the fluid branch results from a growth in the dynami- 
cal correlations. Yet, we show that Xi by itself cannot be 
used as a reliable indirect measure of £4 time-dependence. 

We are grateful to Ludovic Berthier for important dis- 
cussions and insightful comments. 
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FIG. 5. (Color Online) Measurements of X*{t) corrected for 
removal of global effects created while cooling the system. 
Maximum value as well as the time of maximum power law 
diverges (inset) as the termination density p t is approached. 
Maximum value diverges with exponent of 2.0 and the the 
time T4 with exponent of 0.8 

to the dynamics) should be removed. We therefore use 
the measured dynamic correlations Gi(r,t) in order to 
correct the Xi measurements: the exponential decaying 
region of Gi{r, t) is extrapolated to infinite r, and X4 is 
calculated by integrating Gi(r,t) exactly over the short- 
range and adding the contribution of the integral over 
the exponential regime for larger distances. 

The results for the dynamic xa are presented in Figure 
[SJ Peak values of X4 diverges with a power law as pt is ap- 
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